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Abstract 

In  this  paper  we  investigate  the  yield  condition  in  the  mobilization  of  yield- 
stress  materials  in  distensible  tubes.  We  discuss  the  two  possibilities  for 
modeling  the  yield-stress  materials  prior  to  yield:  solid-like  materials  and 
highly-viscous  fluids  and  identify  the  logical  consequences  of  these  two  ap- 
proaches on  the  yield  condition.  As  part  of  this  investigation  we  derive  an 
analytical  expression  for  the  pressure  field  inside  a distensible  tube  with  a 
Newtonian  flow  using  a one-dimensional  Navier-Stokes  flow  model  in  con- 
junction with  a pressure-area  constitutive  relation  based  on  elastic  tube  wall 
characteristics. 

Keywords:  fluid  mechanics;  yield-stress;  yield  condition;  distensible  tube; 
pressure  field. 


1 Introduction 

Many  naturally-occurring  and  synthetic  materials  are  characterized  by  being  yield- 
stress  fluids,  that  is  they  behave  like  solids  under  low  shear  stresses  and  as  fluids  on 
exceeding  a minimum  threshold  yield-stress  [1-5].  Blood  and  crude  oils  are  some 
examples  of  naturally-occurring  yield-stress  fluids  while  many  manufactured  prod- 
ucts used  in  the  food  and  petrochemical  industries;  like  purees,  yoghurts,  slurries, 
drilling  muds,  and  polymeric  solutions;  are  examples  of  synthetic  yield-stress  ma- 
terials [3,  6-20].  As  these  materials  can  reside  and  mobilize  in  rigid  conduits  and 
structures,  such  as  oil  transportation  pipelines  and  geological  porous  formations, 
they  can  also  reside  and  mobilize  in  compliant  conduits  and  structures  like  blood 
vessels  and  living  porous  tissues.  Hence  the  investigation  of  the  circumstances  and 
conditions  under  which  these  materials  yield  and  start  flowing  in  rigid  and  dis- 
tensible containers  is  important  for  modeling  and  analyzing  many  industrial  and 
natural  flow  systems. 
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We  are  not  aware  of  a previous  attempt  to  find  the  condition  for  the  yield  of 
yield-stress  materials  in  distensible  tubes.  Past  investigations  are  generally  focused 
on  the  flow  of  yield-stress  fluids  in  rigid  conduits  and  structures  such  as  ensembles 
of  conduits  and  porous  media  [4,  21-29].  There  are  some  other  studies  (e.g.  [30]) 
that  investigated  the  flow  of  yield-stress  fluids  in  distensible  conduits  but  the  focus 
of  these  studies  is  not  on  the  yield  condition  but  on  the  flow  of  these  materials  as- 
suming they  have  already  reached  a fluid  phase  state  by  satisfying  the  mobilization 
condition. 

To  clarify  the  purpose  of  the  current  investigation,  we  dehne  the  yield  condition 
in  general  terms  that  include  rigid  and  distensible  flow  conduits  as  the  minimum 
pressure  drop  across  the  conduit  that  mobilizes  the  yield-stress  material  and  forces 
it  to  flow  following  an  initial  solid  state  condition.  The  yield  point  is  therefore 
reached  through  increasing  the  pressure  drop  applied  on  the  material  in  its  solid 
state  gradually  or  suddenly.  This  dehnition  is  also  valid  in  principle  for  the  yield 
condition  in  the  bulk  flow  although  this  is  not  of  interest  to  us  in  the  current 
investigation. 

In  the  present  paper  we  investigate  the  yield  condition  in  the  mobilization  of 
yield-stress  materials  through  distensible  tubes  using  an  elastic  model  for  the  tube 
wall  that  correlates  the  pressure  at  a given  axial  coordinate  to  the  corresponding 
cross  sectional  area.  Although  we  assume  a circular  cylindrical  tube  with  a constant 
unstressed  radius  over  its  length,  some  of  the  presented  arguments  at  least  can  be 
generalized  to  include  other  cases  although  we  are  not  intending  to  do  so  in  the 
current  study. 

In  this  investigation,  we  are  only  concerned  with  the  yield  condition  as  it  is, 
without  interest  in  other  issues  related  to  the  flow  of  yield-stress  fluids  in  these 
conduits;  hence  any  other  dynamic  issues  associated  with  the  flow  phase  are  out 
of  scope  of  the  present  paper.  Our  plan  for  the  paper  is  to  introduce  the  one- 
dimensional Navier-Stokes  flow  model  which  is  widely  used  to  describe  the  flow  of 
Newtonian  fluids  in  distensible  tubes.  The  reason  for  using  this  model,  which  is  a 
Newtonian  model  and  not  a yield-stress  model,  is  that  prior  to  yield  the  yield-stress 
materials  according  to  one  approach  behave  as  highly- viscous  Newtonian  fluids. 
We  also  discuss  the  pressure  held  inside  the  tube  under  such  a Navier-Stokes  how 
condition  and  how  this  held  can  be  computed  analytically  and  numerically  because 
it  is  needed  for  identifying  the  yield  condition.  This  will  be  followed  by  a section 
in  which  we  investigate  the  yield  condition  in  detail  considering  two  approaches 
for  modeling  the  yield-stress  materials  prior  to  yield.  A general  assessment  of  our 
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proposed  method  for  identifying  the  yield  condition  with  some  synopsis  conclusions 
will  then  follow. 

2 One- Dimensional  Navier-Stokes  Model 

The  flow  of  Newtonian  fluids  in  a circular  cylindrical  tube  with  length  L and  cross 
sectional  area  A assuming  an  incompressible  laminar  axi-symmetric  flow  with  a 
negligible  gravitational  force  is  described  by  the  following  one-dimensional  Navier- 
Stokes  system  of  mass  and  momentum  conservation  principles 


where  Q is  the  volume  flow  rate,  x is  the  tube  axial  coordinate  along  its  length,  t is 
the  elapsed  time  for  the  flow  process,  p is  the  fluid  mass  density,  a is  the  correction 
factor  for  the  momentum  flux,  p is  the  axial  pressure  along  the  tube  length  which 
is  a function  of  the  axial  coordinate,  and  k is  a viscosity  friction  coefficient  which 
is  normally  defined  hj  k = where  p is  the  ffuid  dynamic  viscosity  [31-34], 

In  this  context  we  assume  a no-slip  at  wall  condition  where  the  velocity  of  the 
fluid  at  the  interface  is  identical  to  the  velocity  of  the  solid  [35-37],  This  Navier- 
Stokes  system  is  normally  supported  by  a constitutive  relation  that  links  the  cross 
sectional  area  at  a certain  axial  location  to  the  corresponding  axial  pressure  in  a 
distensible  tube,  to  close  the  system  in  the  three  variables  Q,  A and  p and  hence 
provide  a complete  mathematical  description  for  such  a flow  in  such  a conduit. 

The  relation  between  the  cross  sectional  area  and  the  axial  pressure  in  a dis- 
tensible tube  can  be  described  by  many  mathematical  models  depending  on  the 
characteristics  of  the  tube  wall  and  its  mechanical  response  to  pressure  such  as 
being  linear  or  non-linear,  and  elastic  or  viscoelastic.  The  following  is  an  example 
of  a commonly  used  pressure-area  constitutive  elastic  relation  that  describes  such 
a dependency 


where  Ao  is  the  reference  cross  sectional  area  corresponding  to  the  reference  pressure 
which  in  this  equation  is  set  to  zero  for  convenience  without  affecting  the  generality 
of  the  results,  A is  the  tube  cross  sectional  area  corresponding  to  the  axial  pressure 


(3) 
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p as  opposite  to  the  reference  pressure,  and  /3  is  the  tube  wall  stiffness  coefficient 
which  is  usually  defined  by 

_ y/nhoE 

where  ho  is  the  tube  wall  thickness  at  the  reference  pressure,  while  E and  ^ are 
respectively  the  Young’s  elastic  modulus  and  Poisson’s  ratio  of  the  tube  wall.  In 
this  context,  we  assume  a constant  ambient  pressure  surrounding  the  tube  that  is 
set  to  zero  and  hence  the  reference  cross  sectional  area  represents  unstressed  state 
with  Ao  being  constant  over  the  axial  coordinate. 

Now,  to  find  the  axial  pressure  field  inside  an  elastic  tube  characterized  by 
the  mechanical  response  of  Equation  3 we  combine  the  pressure-area  constitutive 
relation  with  the  mass  and  momentum  equations  of  Navier-Stokes  system.  In  the 
following  we  assume  a pressure  drop  applied  on  the  tube  by  imposing  two  boundary 
conditions  at  x = 0 and  x = L and  hence  it  is  identified  by  two  pressure  boundary 
conditions,  Pin  and  Pou,  where 


Pin  =p{x  = Q)  > pou  = p{x  = L)  (5) 

Another  clarifying  remark  is  that  for  a sustainable  flow  in  a distensible  tube,  the 
tube  axial  pressure  should  be  a monotonically  decreasing  function  of  its  axial  co- 
ordinate. As  a result,  the  radius  and  the  cross  sectional  area  of  the  tube  are  also 
monotonically  decreasing  functions  of  the  axial  coordinate.  This  fact  is  assumed 
in  most  of  the  following  arguments  although  it  may  not  be  stated  explicitly. 

From  the  pressure-area  constitutive  relation  we  have 

P=^(y^-V^)  (6) 

and  hence 


For  a steady  state  flow  the  time  terms  in  the  Navier-Stokes  system  are  zero  and 
hence  from  the  continuity  equation  (Equation  1)  Q as  a function  of  x is  constant. 
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Consequently,  from  the  momentum  equation  of  the  Navier-Stokes  system  (Equation 
2)  we  obtain 


d f aQ^\  A dp  Q _ 

which  can  be  manipulated  to  obtain  the  axial  pressure  gradient  as  follow 


(9) 


^ dA  V ^ / dp  dx  p dx  A 

aQ'^  dA  dp  A dp  ^ n 
dp  dx^  p dx^  A 


0 


(10) 

(11) 


that  is 


aQ^  ( 2Al  ^ dp  A dp  ^ Q _ r. 


dp  r aQ^  ( 2Al  2A^P\  A 

dx  ^2  I ^2  P + fj  I + p 


(12) 

(13) 


dp 

dx  aQ^  f 2^?  „ _|_  \ A 

[/3^P^  /3  J p 

The  last  expression  can  be  simplihed  to 


(14) 


dp  kq  (yp  + 


(15) 


which  expresses  the  axial  pressure  gradient,  as  a function  of  the  tube  axial 
pressure,  p,  only.  For  a one- dimensional  steady  state  flow,  the  pressure  is  dependent 
on  the  tube  axial  coordinate  only  and  hence  the  partial  derivative  of  Equation  15 
can  be  replaced  with  a total  derivative.  On  separating  the  variables  and  carrying 
out  the  integration,  where  0 < a:  < L and  pou  < P < Pin,  we  obtain 


X = 


2aQ  In  ( ^p  y/Ao 


(3  i^^p  + 


+ C 


(16) 


K bpnQAo 

where  C is  the  constant  of  integration  which  can  be  determined  from  the  boundary 
conditions.  At  a;  = 0 p = pin,  and  hence 
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Pin  “1“ 


c 
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(17) 


The  analytical  solution  of  the  pressure  held,  as  given  by  Equation  16,  dehnes  the 
axial  pressure  held  p{x)  as  an  implicit  function  of  the  axial  coordinate  x.  Obtaining 
the  pressure  held  then  requires  either  the  employment  of  a simple  numerical  solver, 
based  for  example  on  a bisection  method,  or  changing  the  role  of  the  independent 
and  dependent  variables  and  hence  obtaining  x as  a function  of  p where  the  value  of 
pressure  is  constrained  by  the  two  boundary  conditions.  The  latter  scheme,  which 
is  more  convenient  to  use,  leads  to  the  same  result  as  the  hrst  scheme  although  the 
dehning  x-p  points  are  determined  by  the  pressure  points  and  hence  obtaining  a 
smoothly  dehned  pressure  held  may  require  computing  more  points  than  is  needed 
from  the  hrst  scheme.  In  addition  to  the  convenience,  the  solution  obtained  from 
the  second  scheme  is  generally  more  accurate  because  it  produces  exact  solutions 
within  the  available  computing  precision  although  this  extra  accuracy  may  be  of 
little  value  in  practical  circumstances. 

The  analytical  solution  of  Equation  16  was  tested  and  validated  by  numerical 
solutions  from  many  typical  how  examples  using  the  lubrication  approximation 
with  a residual-based  Newton-Raphson  solution  scheme.  This  numerical  method 
for  obtaining  the  pressure  held  and  how  rate  is  based  on  discretizing  the  how 
conduit  into  thin  rings  which  are  treated  as  an  ensemble  of  serially  connected  tubes 
over  which  a mass  conserving  characteristic  how  (in  this  case  the  one-dimensional 
Navier-Stokes  how  in  elastic  cylindrical  tubes)  is  sought  by  forming  a system  of 
simultaneous  equations  derived  from  the  boundary  conditions  for  the  boundary 
nodes  and  the  continuity  of  volumetric  how  rate  for  the  internal  nodes.  This 
system  is  then  solved  numerically  by  an  iterative  non-linear  solution  scheme,  such 
as  Newton-Raphson  method,  to  obtain  the  axial  pressure  held  inside  the  conduit, 
dehned  at  the  boundary  and  internal  nodes,  as  well  as  the  mass-conserving  how 
rate  within  a given  error  tolerance.  This  numerical  solving  scheme  for  obtaining 
the  pressure  held  and  how  rate  is  fully  described  in  [38-40] . A sample  of  the  results 
comparing  the  analytical  to  the  numerical  solutions  for  some  typical  examples  with 
a range  of  tube,  how  and  huid  parameters  are  given  in  Figure  1. 
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Figure  1:  Comparing  the  analytical  solution  of  the  axial  pressure  field  as  given  by 
Equation  16  to  the  numerical  solution  from  the  residual-based  lubrication  approx- 
imation for  the  Navier-Stokes  flow  in  elastic  tubes  with  (a)  R = 0.01,  L = 0.1, 
a = |,  = 738.5,  p = 900,  p = 0.01,  pin  = 2000,  pou  = 0,  and  Q = 2.36708  x 10“^, 

(b)  R = 0.05,  L = 0.75,  a = 1.25,  13  = 3556.0,  p = 1300,  p = 0.025,  pin  = 2500, 
Pou  = 500,  and  Q = 7.34585  x (c)  R = 0.003,  L = 0.045,  a = 1.5,  13  = 70.9, 
p = 1100,  /i  = 0.015,  Pin  = 1900,  Pou  = 300,  and  Q = 4.66119  x 10“®,  and  (d) 
R = 0.1,  L = 0.9,  a = 1.3,  [3  = 23336.4,  p = 800,  p = 0.035,  Pin  = 3100,  Pou  = 600, 
and  Q = 2.49237  x 10“^.  All  the  dimensional  quantities  are  in  standard  SI  units. 
In  all  four  sub-figures,  the  vertical  axis  represents  the  tube  axial  pressure  in  pascals 
while  the  horizontal  axis  represents  the  tube  axial  coordinate  in  meters. 
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3 Yield  Condition 


A general  definition  has  been  given  earlier  for  the  yield  condition.  However,  since 
the  pressure  field  inside  a distensible  tube  and  the  resulting  flow  rate  are  dependent 
on  the  actual  value  of  the  two  boundary  conditions  and  not  only  on  their  difference 
(i.e.  pressure  drop)  [41]  we  define  the  yield  condition  for  distensible  tubes  as  the 
minimum  inlet  pressure  for  a given  outlet  pressure  that  mobilizes  the  yield-stress 
material  and  initiates  a measurable  flow.  The  generally  accepted  condition  for  the 
yield  of  a yield-stress  material  in  a tube  is  given  by  equating  the  magnitude  of  the 
wall  shear  stress  |r^|  to  the  yield-stress  of  the  fluid  Tq  where  is  defined  as  the 
ratio  of  the  force  normal  to  the  tube  axis  F±  to  the  area  of  the  luminal  surface 
parallel  to  this  force  Ay. 

Our  proposal  for  the  condition  that  should  be  satisfied  to  reach  the  yield  point 
is  that  yield  occurs  ijf  all  cross  sectional  points  in  the  tube  along  its  length  reached 
their  yield  condition  simultaneously.  This  will  only  occur  if  the  bottleneck  in  the 
tube  reaches  its  yield  condition  where  the  bottleneck  is  the  location  along  the 
tube  axis  with  the  lowest  magnitude  of  wall  shear  stress.  It  is  obvious  that  this 
is  a necessary  and  sufficient  condition  for  yield  to  occur  because  if  the  bottleneck 
reached  its  yield  condition  with  its  wall  shear  stress  just  exceeding  the  yield-stress 
then  all  other  points  along  the  tube  axis  should  also  satisfy  the  yield  condition. 

For  a rigid  tube  with  a constant  radius  along  its  axial  direction,  Tw  is  given  by 

F^_  ttR^AP  rap  , , 

where  R and  L are  the  tube  radius  and  length  respectively  while  AP  is  the  pressure 
drop  across  the  tube. 

For  a tube  with  a variable  radius  along  its  axial  direction,  including  a distensible 
tube  subjected  to  a pressure  gradient  along  its  axis,  the  wall  shear  stress  is  a 
function  of  the  tube  axial  coordinate  x and  hence  the  wall  shear  stress  is  defined 
locally  as  a function  of  x by  using  inhnitesimal  quantities  of  pressure  drop,  and 
length,  5x,  belonging  to  a thin  cross  sectional  ring  over  which  the  change  in  radius 
is  negligible,  that  is 

Fj_  7rR'^6p  R 6p 

which  in  the  limit  becomes 
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(20) 


R dp 
= 2 * 

Therefore  the  condition  for  the  transition  from  the  solid  state  to  fluid  phase, 
corresponding  to  the  yield  and  flow  initiation,  is  dehned  by  the  following  condition 


R{x) 

dp 

1 to  1 — '1  o ^ 

2 

dx 

= Tn 


(21) 


There  are  two  possibilities  for  modeling  the  yield-stress  materials  before  reach- 
ing their  yield  point:  either  they  are  solid-like  materials  or  they  are  highly- viscous 
fluids.  A detailed  discussion  about  this  issue,  among  other  relevant  issues,  is  given 
in  [29].  Our  proposal  for  identifying  the  yield  point  according  to  each  one  of  these 
two  possibilities  is  discussed  in  the  following  two  subsections. 


3.1  Solid-Like  Materials 


According  to  this  approach,  the  yield-stress  materials  before  reaching  their  yield 
point  behave  like  solids.  Hence,  a logical  assumption  about  the  pressure  held  con- 
hguration  inside  the  tube  prior  to  mobilization  is  to  assume  a linear  pressure  drop 
and  hence  a constant  pressure  gradient.  In  this  case  the  bottleneck  occurs  at  the 
outlet  boundary  point  because  the  radius  of  the  tube  is  monotonically  decreasing 
and  hence  the  outlet  point  is  the  one  with  the  smallest  radius  which  implies  that 
the  wall  shear  stress  as  given  by  Equation  21  is  the  lowest  since  ^ is  assumed  con- 
stant along  the  tube  axis.  Therefore  the  yield  condition  according  to  the  solid-like 
modeling  approach  is  given  by 


R{L)  dp 
2 dx 

where  R{L)  is  the  tube  radius  at  the  outlet  boundary  corresponding  to  x 


(22) 

L. 


3.2  Highly-Viscous  Fluids 

According  to  this  approach,  the  yield-stress  materials  before  reaching  their  yield 
point  are  highly-viscous  fluids  and  hence  they  flow  but  with  an  inhnitesimally  small 
flow  rate.  Our  proposal  in  this  case  is  that  the  material  prior  to  yield  should  be 
modeled  as  a Newtonian  fluid  with  a very  high  viscosity.  The  assumption  of  a 
Newtonian  fluid  is  justihed  by  the  fact  that  at  such  regimes  of  very  low  rate  of 
deformation  the  material  is  at  its  low  Newtonian  plateau  since  all  non-Newtonian 
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rheological  effects  are  induced  only  by  a measurable  deformation  [29] . We  therefore 
use  the  one-dimensional  Navier-Stokes  flow  model  to  identify  the  pressure  held  prior 
to  yield  and  hence  the  yield  condition.  This  one-dimensional  Navier-Stokes  how 
model  is  fully  described  in  [32,  33].  It  is  also  outlined  in  section  2 with  a derivation 
of  an  analytical  expression  for  the  axial  pressure  as  an  implicit  function  of  the  tube 
axial  coordinate.  The  determination  of  the  axial  pressure  held  is  obviously  needed 
for  the  identihcation  of  the  yield  condition  because  both  R and  ^ are  dependent 
on  the  pressure  held. 

Our  general  approach  for  identifying  the  yield  condition  assuming  a highly- 
viscous  huid  prior  to  yield  is  to  test  the  axial  pressure  held  as  determined  by  the 
given  inlet  and  outlet  boundary  conditions.  The  purpose  of  this  test  is  to  identify 
the  bottleneck  by  determining  the  lowest  wall  shear  stress  in  magnitude  along  the 
tube  and  hence  verifying  if  the  bottleneck  has  reached  its  minimum  threshold  which 
equals  the  yield-stress  value,  Tq.  As  seen  in  Equation  21,  |r^|  is  dependent  on  the 
pressure  gradient  and  the  tube  radius  as  functions  of  the  axial  coordinate.  Hence 
a knowledge  of  these  two  quantities  is  required. 

There  are  two  ways  for  obtaining  the  pressure  gradient  and  the  tube  radius. 
The  hrst  way  is  to  solve  the  pressure  held  numerically  using  the  aforementioned 
residual-based  lubrication  approach.  The  axial  pressure  gradient  is  then  calculated 
from  the  numerical  solution  of  the  pressure  held  either  analytically  from  Equation 
15  or  numerically  using  for  instance  a hnite  diherence  method.  On  knowing  the 
axial  pressure,  the  tube  radius  as  a function  of  the  axial  coordinate  can  also  be 
computed  from  the  pressure-area  constitutive  relation  (Equation  7).  The  second 
way  is  to  use  the  analytical  solution  of  the  pressure  held,  which  is  derived  in  section 
2,  and  hence  evaluate  the  axial  pressure  gradient  and  axial  radius  as  in  the  case 
of  numerical  solution  of  the  pressure  held,  although  the  analytical  evaluation  of 
the  pressure  gradient  from  Equation  15  may  not  be  possible  unless  Q is  obtained 
from  another  method.  Both  of  these  ways  practically  produce  the  same  result 
if  a sufficiently  hne  mesh  is  used  for  discretizing  the  tube  in  the  residual-based 
lubrication  technique. 

As  indicated  early,  in  the  case  of  using  the  analytical  expression  of  Equation  16 
a simple  numeric  solver,  based  for  instance  on  a bisection  method,  can  be  used  to 
evaluate  the  axial  pressure  which  is  given  as  an  implicit  function  of  the  tube  axial 
coordinate.  An  alternative  and  more  convenient  way  that  does  not  require  the 
employment  of  a numerical  solver  is  to  use  Equation  16  to  hnd  the  axial  coordinate 
a:  as  a function  of  the  axial  pressure  p,  which  is  as  good  for  evaluating  the  axial 
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pressure  field  as  the  other  way  round.  The  use  of  one  of  these  methods  or  the  other 
(i.e.  numerical  or  analytical  with  or  without  the  employment  of  a numerical  solver) 
is  a matter  of  choice  and  convenience.  All  our  numerical  experiments  produced 
essentially  the  same  results  on  using  these  different  methods  and  techniques. 

However,  it  can  be  shown  that  the  bottleneck,  according  to  the  highly-viscous 
fluid  approach,  is  always  at  the  inlet  boundary  and  hence  the  use  of  the  numerical 
and  analytical  methods  is  unnecessary  apart  from  providing  the  one-sided  pressure 
gradient  at  the  inlet  since  the  inlet  radius  is  already  known  from  the  inlet  boundary 
condition.  The  justihcation  of  this  is  explained  in  the  following  paragraphs. 

As  stated  before,  for  a steady  state  flow  the  time  terms  in  the  Navier-Stokes 
system  are  zero  and  hence  from  the  continuity  equation  (Equation  1)  Q as  a func- 
tion of  X is  constant.  Hence,  from  the  momentum  equation  of  the  Navier-Stokes 
system  (Equation  2)  we  obtain 


d ( aQ^\  ,^dp  Q _ 


p 

r 

A 

(aQ^\\ 

A 

L A ^ 

dx 
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Now  since  A is  a monotonically  decreasing  function  of  x 
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In  the  following  we  will  show  that  / = |r^^,|  is  a monotonically  increasing  function 
of  X by  demonstrating  that  the  hrst  derivative  of  / is  positive,  i.e.  ^ > 0. 
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Now,  is  either  negative  or  positive.  If  it  is  negative  then  the  condition  given 
by  Eqnation  42  is  always  satisfied  because  the  right  hand  side  is  strictly  positive. 
However,  if  it  is  positive  then  since 


B 

C 


K 


2aQ 


(43) 


we  can  show  that  in  all  practical  situations  the  first  term  on  the  right  hand  side 
is  greater  than  the  left  hand  side  regardless  of  the  magnitude  of  the  second  term 
on  the  right  hand  side.  The  reason  is  that  prior  to  yield  Q is  infinitesimally  small, 
and  therefore  ^ is  very  large  in  magnitude  and  positive  in  sign.  Hence  for  any 
tangible  pressure  drop  that  makes  |^|  hnite,  which  is  automatically  satisfied  for 
any  tangible  yield-stress,  the  first  term  on  the  right  hand  side  is  very  large.  There- 
fore, for  all  practical  purposes  where  the  curvature  of  the  tube,  as  quantified  by 
the  second  derivative  of  R with  respect  to  x,  is  constrained  within  the  physical  lim- 
its the  condition  given  by  Equation  42  will  be  satisfied.  This  means  that  in  both 
cases  the  hrst  derivative  of  / is  greater  than  zero  and  hence  the  absolute  value  of 
the  wall  shear  stress,  is  a monotonically  increasing  function  of  the  tube  axial 
coordinate  x. 

As  a result,  the  bottleneck  according  to  the  highly- viscous  fluid  approach  will 
always  be  at  the  inlet  boundary.  This  Ending  is  confirmed  by  numerical  experiments 
where  all  the  results  indicate  that  the  bottleneck  is  at  the  inlet  boundary.  It  should 
be  remarked  that  it  can  be  demonstrated  that  for  the  one-dimensional  Navier- 
Stokes  flow  in  an  elastic  tube  whose  mechanical  response  is  described  by  Equation 
3 it  is  always  the  case  that  < 0,  i.e.  the  axial  dependency  of  the  radius,  area 
and  pressure  concave  downward  (refer  for  instance  to  Figure  1).  However,  we  will 
not  prove  this  here  to  avoid  unnecessary  details  because,  as  demonstrated  already, 
this  is  not  needed  for  establishing  our  case. 

To  sum  up,  the  identification  of  the  bottleneck  through  the  inspection  of  the 
axial  pressure  field  by  using  the  numerical  or  the  analytical  solution  is  the  generally 
valid  method  and  hence  it  should  provide  a definite  and  reliable  answer  for  deter- 


13 


mining  the  yield  condition.  However,  as  demonstrated  in  the  previous  paragraphs, 
at  least  for  all  the  practical  purposes  in  which  the  distensible  tube  is  not  subjected 
to  extremely  large  pressure  boundary  conditions  which  invalidate  the  pressure-area 
relation,  the  bottleneck  and  hence  the  yield  condition  can  be  identihed  from  just 
inspecting  the  magnitude  of  the  wall  shear  stress  at  the  inlet  boundary,  which  as  ex- 
plained earlier  only  requires  computing  the  one-sided  pressure  gradient  numerically 
or  analytically.  In  all  cases,  tests  can  be  carried  out  to  verify  that  the  bottleneck 
is  actually  at  the  inlet  boundary. 

4 General  Assessment 

The  method  proposed  in  this  paper  for  the  identihcation  of  the  yield  condition  in 
the  mobilization  of  yield-stress  materials  in  distensible  tubes  is  based  on  a number 
of  simplifying  assumptions  regarding  the  fluid,  flow  and  tube.  The  purpose  of  this 
investigation  is  to  derive  the  yield  condition  from  logical  reasoning  through  the 
employment  of  the  principles  of  mechanics,  rheology  and  fluid  dynamics.  Many 
practical  factors  that  usually  arise  in  real  life  can  change  this  condition  and  hence 
require  a deeper  inspection  to  identify  the  yield  condition  more  rigorously.  The  real 
yield-stress  flow  systems  are  more  physically  complex  than  the  description  provided 
by  our  model  or  any  other  model  since  all  these  models  are  based  on  mathematical 
idealizations  and  physical  simplihcations.  However,  we  believe  this  investigation 
can  provide  the  ground  for  further  investigations  in  the  future  through  which  more 
complex  factors  can  be  incorporated  in  the  modeled  yield-stress  flow  systems  to 
provide  better  predictions  for  the  yield  condition. 

5 Conclusions 

The  investigation  in  the  present  paper  lead  to  the  proposal  of  a method  for  iden- 
tifying the  yield  condition  in  the  mobilization  of  yield-stress  materials  through 
deformable  cylindrical  tubes  with  a constant  unstressed  cross  sectional  area  along 
their  axial  length.  Two  possible  scenarios  for  modeling  the  yield-stress  materials 
prior  to  yield  have  been  considered:  a solid-like  approach  and  a highly-viscous  fluid 
approach.  The  logical  consequences  of  these  two  approaches  were  derived  using  the 
principles  of  rheology  and  fluid  dynamics  as  summarized  in  the  rheological  char- 
acteristics of  the  deformed  material,  Navier-Stokes  system  and  the  distensibility 
characteristics  of  the  tube  wall.  General  mechanical  principles  were  also  employed 
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in  this  investigation  with  regard  to  the  solid-like  approach. 

Although  the  derived  yield  condition  is  based  on  several  simplifying  assumptions 
with  regard  to  the  type  of  flow,  fluid  and  tube  characteristics  and  hence  it  applies  to 
a rather  idealized  yield-stress  flow  system,  the  proposed  model  can  serve  as  a first 
step  for  more  elaborate  yield-stress  models  that  incorporate  more  physical  factors 
that  influence  the  yield  condition. 

In  the  course  of  this  investigation,  an  analytical  expression  linking  the  axial 
pressure  field  inside  a compliant  tube  having  elastic  properties  to  its  axial  coor- 
dinate along  its  length  has  also  been  derived  and  validated  by  a residual-based 
numerical  method  from  the  lubrication  theory. 
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Nomenclature 


a correction  factor  for  axial  momentum  flux 
13  stiffness  coefficient  in  pressure-area  relation 

K viscosity  friction  coefficient 

fj,  ffuid  dynamic  viscosity 

p ffuid  mass  density 

^ Poisson’s  ratio  of  tube  wall 

To  yield-stress 

Tw  shear  stress  at  tube  wall 

A tube  cross  sectional  area 

Ain  tube  cross  sectional  area  at  inlet 

Ao  tube  cross  sectional  area  at  reference  pressure 

Aou  tube  cross  sectional  area  at  outlet 

A\\  luminal  area  parallel  to  tube  axis 

E Young’s  elastic  modulus  of  tube  wall 

F_[_  normal  force 

ho  tube  wall  thickness  at  reference  pressure 
L tube  length 

p pressure 

Pin  inlet  boundary  pressure 

Pou  outlet  boundary  pressure 

5p  inhnitesimal  pressure  drop 

AP  pressure  drop  across  the  tube 
Q volumetric  flow  rate 

R tube  radius 

t time 

X tube  axial  coordinate 

5x  inhnitesimal  length  along  tube  axis 
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